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Abstract 

Recent deep inelastic scattering and Drell-Yan experiments have revealed a surpris- 
ingly large asymmetry between the up and down sea quark distributions in the nucleon. 
This result strongly suggests that the mesonic degrees of freedom play an important 
role in the description of the parton distributions in hadrons. In this article, we review 
the current status of the flavor asymmetry of the nucleon sea. The implications of the 
pionic models as well as possible future measurements are also discussed. 
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1. INTRODUCTION 

The first direct evidence for point-like constituents in the nucleons came from the ob- 
servation of scaling phenomenon in Deep-Inelastic Scattering (DIS) experiments ]l[ |2| at 
SLAC. These point-like charged constituents, called partons, were found to be spin- 1/2 
fermions. These partons were initially identified as the quarks in the Constituent Quark 
Models (CQM). However, it was soon realized that valence quarks alone could not account 
for the large enhancement of the cross sections at small Bjorken-x, the fraction of nucleon 
momentum carried by the partons. These swarm of low-momentum partons, dubbed wee- 
partons by Feynman ||, was interpreted as the quark and antiquark sea of the nucleon. 
DIS experiments therefore provided the first evidence for the existence of antiquarks in the 
nucleon. 

The observation of partons in DIS experiments paved the road to the formulation of 
Quantum Chromodynamics (QCD) as the theory for strong interaction. Nevertheless, the 
exact form of the parton distribution functions can not be deduced from perturbative QCD. 
Like many other static properties of hadrons, the parton distribution functions belong to 
the domain of non-perturbative QCD. In spite of great progress |4j] made in Lattice Gauge 
Theory (LGT) in treating the bound-state properties of hadrons, it remains a challenge to 
predict the parton distributions using LGT. 

Until parton distributions can be readily calculated from first principles, they are best 
determined from experiments. Electroweak processes such as DIS and lepton-pair production 
provide the cleanest means to extract information on the parton distributions. There are at 
least two reasons why it is important to measure the parton distribution functions. First, 
the description of hard processes in high energy interactions requires parton distribution 
functions as an essential input. Second, many aspects of the parton distributions, such as 
sum rules, scaling- violation, asymptotic behaviors at large and small x, and flavor and spin 
structures, can be compared with the predictions of perturbative as well as non-perturbative 
QCD. 

In this article, we review the status of our current knowledge on the flavor dependence of 
the sea quark distributions in hadrons. The recent observation of a striking flavor asymmetry 
of the nucleon sea has profound implications on the importance of meson degrees of freedom 
for the description of parton substructures in the nucleon. In Section 2, we review the early 
studies of the nucleon sea in DIS and lepton-pair production. The crucial recent experiments 
establishing the up/down flavor asymmetry of the nucleon sea are discussed in Section 3. 
The various theoretical models for explaining the u/d asymmetry are described in Section 4. 
The implications of meson-cloud models on other aspects of the parton structure functions 
are discussed in Section 5. Finally, we present future prospects in Section 6, followed by 
conclusion in Section 7. 

2. EARLY STUDIES OF THE NUCLEON SEA 
2.1 Deep Inelastic Scattering 

Although scaling behavior in inelastic electron scattering was predicted by Bjorken || 
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based on the framework of current algebra, its confirmation by the SLAC experiments still 
came as a major surprise. A simple and intuitive picture for explaining the scaling behavior 
is the parton model advanced by Feynman |3|, ||. In this model, the electron- nucleon deep- 
inelastic scattering is described as an incoherent sum of elastic electron-parton scattering. 
However, the nature of the partons within the nucleon was not specified by Feynman. Some 
authors [0, || |9], [T(| speculated that the partons are the 'bare nucleon' plus the pion cloud, 



while others [JTT] , |12| , |13| believed they are the quarks introduced by Gell-Mann [L4j and 
Zweig ||15|| . The latter scenario was strongly supported by the measurement of R, ratio of 
the longitudinally over transversely polarized photon cross sections, showing the spin- 1/2 
character of the partons. 

Evidence for quark-antiquark sea in the nucleon came from the observation that the 
structure function F 2 (x) approaches a constant value as x — > JTB|. If the proton is made up 
of only three quarks, or any finite number of quarks, F 2 (x) is expected to vanish as x — > 0. 
Bjorken and Paschos [|l l| therefore assumed that the nucleon consists of three quarks in 
a background of an indefinite number of quark-antiquark pairs. Kuti and Weisskopf |T2| 
further included gluons among the constituents of nucleons in order to account for the missing 
momentum not carried by the quarks and antiquarks alone. 

The importance of the quark-antiquark pairs in the nucleon is in sharp contrast to the 
situation for the atomic system, where particle-antiparticle pairs play a relatively minor role 
(such as the polarization of the vacuum). In strong interactions, pairs of virtual fermions or 
virtual Goldstone bosons are readily produced as a result of the relatively large magnitude 
of the coupling constant a s , and they form an integral part of the nucleon's structure. 

Neutrino-induced DIS experiments allowed the separation of sea quarks from the valence 
quarks. Recall that 

F 2 p (x) = 2xJ2 [Qi( x ) + Qi( x )] > 



F» p (x) = 2j2Mx)-q t (x)] = 2j2q. 



vl x) 



where i denotes the flavor of the quarks. Note that the valence quark distribution is defined 
as the difference of the quark and antiquark distributions, qf(x) = q%(x) — q~i(x). Eq. |T] shows 
that the valence quark distribution is simply F^ p (x)/2, while the sea quark distribution is 
given by F 2 p (x)/2x — F^ p (x)/2. The F 2 (x) and F 3 (x) data from the CDHS experiment ]T7J 
clearly showed that the valence quark distributions dominate at x > 0.2, while the sea quarks 
are at small x. 

The earliest parton models assumed that the proton sea was flavor symmetric, even 
though the valence quark distributions are clearly flavor asymmetric. Inherent in this as- 
sumption is that the content of the sea is independent of the valence quark's composition. 
Therefore, the proton and neutron are expected to have identical sea-quark distributions. 
The flavor symmetry assumption was not based on any known physics, and it remained to 
be tested by experiments. Neutrino- induced charm production experiments ]TJ, [L9|, which 
are sensitive to the s — > c process, provided strong evidences that the strange-quark content 
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of the nucleon is only about half of the up or down sea quarks. Such flavor asymmetry is 
attributed to the much heavier strange-quark mass compared to the up and down quarks. 
The similar masses for the up and down quarks suggest that the nucleon sea should be nearly 
up-down symmetric. A very interesting method to check this assumption is to measure the 
Gottfried integral in DIS, as discussed next. 



2.2 Gottfried Sum Rule 

In 1967, Gottfried studied electron-proton scattering assuming that the proton consists of 
three constituent quarks |20] . He showed that the total electron-proton cross section (elastic 



plus inelastic) is identical to the Rutherford scattering from a point charge. Gottfried derived 
a sum rule 

if = r FZ{x, Q 2 )/x dx = Y.(Q P if = 1. (2) 
Jo T 

where Q\ is the charge of the zth quark in the proton. Gottfried expressed great skepticism 
that this sum rule would be confirmed by the forthcoming SLAC experiment by stating "I 
think Prof. Bjorken and I constructed the sum rules in the hope of destroying the quark 
model" . Indeed, Eq. || was not confirmed by the experiments, not because of the failure 
of the quark model, but because of the presence of quark-antiquark sea. In fact, the infinite 
number of the sea partons makes If diverge. A closely related sum rule, now called the 
Gottfried Sum Rule (GSR), avoids this problem by considering the difference of the proton 
and neutron cross sections, namely, 

if - i 2 n = [\f 2 p (x, Q 2 ) - F*{z, Q 2 )]/x dx = Y.m? ~ (Q- f] = 1/3- (3) 
Jo T 

In deriving Eq. [| it was assumed that the sea quarks in the proton and neutron are related 
by charge symmetry (see Section 4.5). 

Soon after the discovery of scaling in electron-proton DIS, electron-deuterium scattering 
experiments were carried out to extract the electron-neutron cross sections. The comparison 
of e — p with e — n data was very important for distinguishing early competing theoretical 
models |16[]. These data also allowed a first evaluation |22| of the Gottfried integral in 1970. 



The first result for the Gottfried integral was 0.19, considerably less than 1/3. Note that 
the data only covered x > 0.08, and it was assumed that Ff — F 2 ™ follows Regge behavior 
(proportional to x 1 ^ 2 ) in the unmeasured small- x region. Due to the 1/x factor in the 
integrand, the small-x region could have potentially large contributions to the Gottfried 
integral. Moreover, it was not clear if F% — F£ would indeed follow the Regge behavior at 
small x, and if so, at what value of x. By 1973, new data were available down to x — 0.05 
and the Gottfried integral was evaluated to be 0.28 [p^] , considerably larger than the first 
result. It should be pointed out that these data were taken at relatively low values of Q 2 . 
Furthermore, Q 2 varies as a function of x. 

Although the large systematic errors associated with the unmeasured small- x region 
prevented a sensitive test of the GSR, Field and Feynman Pi nevertheless interpreted the 
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early SLAC data as strong indications that GSR is violated and that the u and d distributions 
in the proton are different. The relationship between the Gottfried integral and the d/u 
asymmetry is clearly seen in the parton model, namely, 

I [FUx,Q 2 )-F^(x,Q 2 )]/x dx = - + - ( l [u{x)-d(x)]dx. (4) 
Jo 3 3 Jo 

Eq. [| clearly shows that the early SLAC data on GSR implied d > u, at least for certain 
region of x. Field and Feynman further suggested that Pauli blocking from the valence 
quarks would inhibit the uu sea more than the dd sea, hence creating an asymmetric nucleon 
sea. 

The SLAC DIS experiments were followed by several muon-induced DIS experiments at 
Fermilab and at CERN. Using high energy muon beams, these experiments reached much 
larger values of Q 2 and they clearly established |2l| 27) the scaling-violation phenomenon 
in DIS. The Gottfried integral was also evaluated in muon DIS experiments |2S], |29|. Figure [I] 
compares the data from the European Muon Collaboration (EMC) |28| with earlier electron 



data from SLAC [ 3(J . The coverages in x are similar in these two experiments, even though 
the Q 2 values covered by the EMC are much larger. Figure [I] shows that F% — F£ from 
EMC tend to shift towards smaller x relative to the SLAC data, in qualitative agreement 
with the QCD Q 2 -evolution. The Gottfried integral determined from the EMC experiment 
is 0.235 + 0.110 - 0.099, consistent with that from the SLAC, but still lower than 1/3. 

Despite the fact that all measurements of Gottfried integral consistently showed a value 
lower than 1/3, the large systematic errors prevented a definitive conclusion. As a result, 
all parametrizations |JT], |3B], |55|] of the parton distributions based on global fits to 



existing data before 1990 assumed a symmetric u, d sea. As discussed later, compelling 
evidence for an asymmetric light-quark sea awaited new experimental inputs. 



2.3 Drell-Yan Process 

The first high-mass dilepton production experiment [p6] was carried out at the AGS in 



1969, soon after scaling was discovered at SLAC. Drell and Yan |37|] interpreted the data 
within the parton model, in which a quark-antiquark pair annihilate into a virtual photon 
decaying subsequently into a lepton pair. This simple model was capable of explaining 
several pertinent features of the data, including the overall magnitude of the cross sections, 
the scaling behavior of the cross sections, and the polarization of the virtual photons. The 
high-mass continuum lepton-pair production is therefore called the Drell-Yan (DY) process. 

Since the underlying mechanism for the DY process involves the annihilation of a quark 
with an antiquark, it is not surprising that this process can be used to probe the antiquark 
contents of the beam or target hadrons. In the parton model, the DY cross section is given 
by 

d^dx~ F = IZs ( Xl I x 2 ) ? ZfoMf'M + (5) 

Here q a (x) are the quark or antiquark structure functions of the two colliding hadrons eval- 
uated at momentum fractions x\ and x%- The sum is over quark flavors. In addition, one 
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has the kinematic relations 



r = X1X2 = M 2 /s, 

xp = X\ — x 2 , (6) 

where M is the invariant mass of the lepton pair and s is the square of the center-of-mass 
energy. The cross section is proportional to a 2 , indicating its electromagnetic character. 

Eq. [5] shows that the antiquark distribution enters as a multiplicative term in the DY 
cross section rather than an additive term in the DIS cross section. Hence, the antiquark 
distributions can be sensitively determined in the DY experiments. The dimuon data from 
the FNAL E288, in which the T resonances were discovered |38| , were analysed |39| to 
extract the sea quark distribution in the nucleon. By assuming a flavor-symmetric nucleon 
sea, namely, u(x) = d(x) = s(x) = Sea(x), the dimuon mass distribution obtained in 400 
GeV proton-nucleus interaction was described by xSea(x) = 0.6a; -10 . In a later analysis (|0| 
containing additional data at 200 and 300 GeV, E288 collaboration found that a much better 
fit could be obtained with an asymmetric sea, namely, 

u(x) = (1 - x) 3A8 d(x), s(x) = (u(x) + d(x))/A. (7) 



The need for an asymmetric u and d was also manifested in the E288 da/dy data ||0| at 
y — 0, where y is the center-of-mass rapidity. For p + A collision, da/dy at y — is expected 
to be positive due to the excess of u over d valence quarks in the proton. The E288 data 
showed that the slopes are indeed positive, but larger than expected from a flavor symmetric 
sea. A surplus of d over u in the proton sea would lead to more positive slope in agreement 
with the data [f|(| . 



The FNAL E439 collaboration |4I[] studied high mass dimuons produced in p + W in- 
teraction at 400 GeV. Their spectrometer covered a considerably larger range in xp than 
E288. They again found that an asymmetric sea, u[x) = (1 — x) 2,5 d(x), can well describe 
their data. 

With all the tantalizing evidence for an asymmetric sea from DIS and DY experiments, 



it is curious that all global analyses pfl, 32, 33, p4l Rf| of parton distributions in the 1980's 



still assumed a symmetric light-quark sea. This probably reflected the reluctance to adopt 
an unconventional description of the nucleon sea without compelling and unambiguous ex- 
perimental evidences. As discussed in the next Section, such evidence became available in 
the 1990's. 



3. RECENT EXPERIMENTAL DEVELOPMENTS 
3.1 NMC Measurements of the Gottfried Integral 



After the discovery of the so-called 'EMC effect' which showed that the parton 
distribution in a nucleus is different from that in the deuteron, the EMC detector system 
was modified by the New Muon Collaboration (NMC) to study in detail the EMC effect. 
Special emphases were placed by the NMC on the capability to explore the low- a; region 



s 



where the 'shadowing effect' is important, and on accurate measurement of cross section 
ratios |44 |. By implementing a 'small angle' trigger which extended the scattering angle 
coverage down to 5 mrad, the lowest value of x reached by NMC was ~ 0.001. The NMC 
also placed two targets in the muon beam, allowing DIS data from two different targets to be 
recorded simultaneously, thus greatly reducing the beam flux normalization uncertainty. To 
account for the different geometric acceptances for events originating from the two targets, 
half of the data were taken using a different target configuration where the locations of the 
two targets were interchanged. 



The upgraded NMC detectors allowed a definitive study of the shadowing effect at low 
x. Moreover, they enabled a much more accurate determination of the Gottfried integral. 
Figure |2] shows the Ff — F£ reported by the NMC in 1991 |]46| , in which the smallest x 
reached (0.004) was significantly lower than in previous experiments. Taking advantage of 



their accurate measurements of the F 2 n / F : 



p ratios, NMC used the following expression to 



evaluate F 2 — Fg, namely, 



Ff - F? 



2 F? 



1 -*?/*?)/(! + 



The ratio F^/Ff = 2F 2 d /Ff — 1 was determined from NMC's F^/Ff measurement, while 



F 2 was taken from a fit to previous DIS experiments at SLAC [J47| , Fermilab [27|, and 



CERN 



49 



The value of the Gottfried integral for the measured region at Q 2 = 4 GeV 2 
is 5 G (0.004-0.8) = 0.227 ± 0.007 (stat) ±0.0U(syst). The contribution to S G from x > 0.8 
was estimated to be 0.002 ± 0.001. Assuming that F 2 — FT at x < 0.004 behaves as ax b , 
NMC obtained Sg(0 — 0.004) = 0.011 ± 0.003. Summing the contributions from all regions 
of x, NMC obtained Sq = 0.240 ±0.016, which was significantly below 1/3. This represented 
the best evidence thus far that the GSR was indeed violated. 



In 1994, NMC reevaluated [HI the Gottfried integral using a new F^ parametrization 



from Ref. [51| and newly determined values of F^/Ff. The new F^ parametrization include 
g7|, BCDMS [(H, as well as NMC's own measurement [[51]]. The new 
F 2 ™ are shown in Figure |2|. Note that the 1994 values are slightly larger 
(smaller) than the 1991 values at small (large) x. The new evaluation gave £g(0. 004 — 0.8) = 
0.221 ± 0.008(stat) ± 0M9(syst), and the Gottfried integral became 



data from SLAC 
NMC results for Ff 



S G = 0.235 ± 0.026. 



(9) 



This value is consistent with the earlier number. The new systematic error is larger than the 
old one, reducing somewhat the overall significance of the NMC measurement. Nevertheless, 
the violation of the GSR is established at a 4<r level. 



More recently, NMC published their final analysis of F^/Ff |2| and F 2 d 0, |54|. This 
analysis includes the 90 and 280 GeV data taken in 1986 and 1987, as well as the 1989 data 
at 120, 200 and 280 GeV. The 1989 data were not used in the earlier evaluations |y| 0] of 
the Gottfried integral. Based on these new data, NMS reported a Gottfried integral in the 
interval 0.004 < x < 0.8 of 0.2281 ± 0.0065(stat) at Q 2 = 4 GeV 2 |52j. This agrees within 
statistical errors with previous NMC results MR, |50| . 
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QCD corrections to various parton-model sum rules have been reviewed by Hinchliffe 
and Kwiatkowski |55fl. The a, and a 2 corrections to the Gottfried sum have been calcu- 



lated 0, |5^] and found to be very small (roughly 0.4 % each at Q 2 = 4 GeV 2 ). Therefore, 
QCD corrections can not account for the large violation of GSR found by the NMC. Al- 
though perturbative QCD predicts a weak Q 2 dependence for the Gottfried sum, it has been 
suggested |58], |59| that due to the non-perturbative origin of the d, u asymmetry the Q 2 
dependence of the Gottfried sum will be anomalous between 1 and 5 GeV 2 . This interesting 
suggestion remains to be tested by DIS experiments. 

3.2 E772 Drell-Yan Experiment 

The main physics goal of the Fermilab experiment 772 was to examine the origin of the 



EMC effect. Among the many theoretical models which explain the EMC effect ||60|| , the pion- 
excess model |HJ predicted a large enhancement of antiquark content due to the presence of 
additional meson cloud in heavy nuclei. This prediction could be readily tested by measuring 
the nuclear dependence of proton-induced DY cross sections. Using an 800 GeV proton beam, 
the DY cross section ratios of a number of nuclear targets (C, Ca, Fe, W) were measured |)2 



relative to deuterium with high statistics over the region 0.04 < x 2 < 0.35. While the 
enhancement of antiquark contents predicted by the pion-excess model was not observed, 
the E772 results are in good agreement with the prediction of the rescaling model [[kJ . 

Information on the d/u asymmetry has also been extracted from the E772 data []64 |. At 



Xf > 0.1, the dominant contribution to the proton-induced DY cross section comes from the 
annihilation of u quark in the projectile with the u quark in the target nucleus. It follows 
that the DY cross section ratio of a non-isoscalar target (like 1 H, Fe, W) over an isoscalar 
target (like 2 H, 12 C, m Ca) is given as 

RaM = v A (x 2 )/a IS (x 2 ) « 1 + [(N-Z)/A}[(l-u(x 2 )/d(x 2 ))/(l + u(x 2 )/d(x 2 ))}, (10) 

where x 2 is the Bjorken-x of the target partons, IS stands for isoscalar, and N, Z and A refer 
to the non-isoscalar target. The (N — Z) / A factor in Eq. [K] shows that the largest sensitivity 
to the d/u could be obtained with a measurement of cr(p + p)/cr(p + d). Nevertheless, for W 
target (N — Z)/A = 0.183 and the E772 aw/vis data could still be used to study the d, u 
asymmetry. 

Figure |3| shows the E772 DY cross section ratios from the neutron-rich W target over 
the isoscalar targets, 2 H and l2 C . Corrections have been applied to the two data points 
at x 2 < 0.1 to account for the nuclear shadowing effects in C and W []64| . The E772 



data provided some useful constraints on the d/u asymmetry. In particular, some early 
parametrizations of large d/u asymmetry were ruled out. Despite the relatively 



large error bars, Figure ^| shows R > 1.0 at x 2 > 0.15, suggesting that d > u in this region. 

E772 collaboration also presented the DY differential cross sections for p + d at (M w ) = 
8.15 GeV. As shown in Figure ^, the DY cross sections near xp = are sensitive to d/u 
and the data disfavor a large d/u asymmetry. Figure f| also shows that a recent parton 
distribution function, MRST |67[| , which has modest d/u asymmetry, is capable of describing 
the p + d differential cross sections well (see Section 3.6). 
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While the E772 data provide some useful constraints on the values of d/u, it is clear 
that a measurement of <tdy(p + d)/aDY{p + p) is highly desirable. The (N — Z)/A factor 
is now equal to —1, indicating a large improvement in the sensitivity to d/u. Moreover, 
the uncertainty arising from nuclear effects would be much reduced. It has been pointed 



out |68fl that d/u asymmetry in the nucleon could be significantly modified in heavy nu- 
clei through recombination of partons from different nucleons. Clearly, an ideal approach 
is to first determine d/u in the nucleon before extracting the d/u information in heavy nuclei. 

3.3 NA51 Drell-Yan Experiment 

Following the suggestion of Ellis and Stirling [B5|, the NA51 collaboration at CERN 



carried out the first dedicated dimuon production experiment to study the flavor structure 



of the nucleon sea [ S9| . Using a 450 GeV/c proton beam, roughly 2800 and 3000 dimuon 
events with > 4.3 GeV have been recorded, respectively, forp + p and p + d interaction. 
The spectrometer setting covers the kinematic region near y = 0. At y — 0, the asymmetry 
parameter, Aoy, is given as 

_ gPP - gPn ^ (4V _ 1)(As _ 1) + (Ay _ 1)(Aj _ 1) 

DY ~ aPP + aP- ~ (4Ay + l)(A s + l) + (Ay + l)(A s + l)' 1 ' 



where Ay = uy/dy and X s = u/d. In deriving Eq. [TTJ, the negligible sea-sea annihilation 
was ignored and the validity of charge symmetry was assumed. At x = 0.18, Ay ~ 2 and 
according to Eq. [TT] Aoy = 0.09 for a symmetric sea, A s = 1. For an asymmetric d > u sea, 
Ady would be less than 0.09. 

From the DY cross section ratio, a pp /a pd , NA51 obtained A DY = -0.09 ± 0.02(stat) ± 
0.025(syst). This then led to a determination of u/d = 0.51 ±0.04(stat) ± (syst) at x = 0.18 
and (M w ) = 5.22 GeV. This important result established the asymmetry of the quark sea 
at a single value of x. What remained to be done was to map out the x-dependence of 
this asymmetry. This was subsequently carried out by the Fermilab E866 collaboration, as 
described in the next Section. 



3.4 E866 Drell-Yan Experiment 

At Fermilab, a DY experiment (E866) aimed at a higher statistical accuracy with a 
much wider kinematic coverage than the NA51 experiment was recently completed |[70| . 
This experiment measured the DY muon pairs from 800 GeV proton interacting with liquid 
deuterium and hydrogen targets. A proton beam with up to 2 x 10 12 protons per 20 s spill 
bombarded one of three identical 50.8 cm long cylindrical target flasks containing either 
liquid hydrogen, liquid deuterium or vacuum. The targets alternated every few beam spills 
in order to minimize time- dependent systematic effects. The dimuons accepted by a 3-dipole 
magnet spectrometer were detected by four tracking stations. An integrated flux of 1.3 x 10 17 
protons was delivered for this measurement. 

Over 330,000 DY events were recorded in E866, using three different spectrometer settings 
which covered the regions of low, intermediate and high mass muon pairs. The data presented 
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here are from the high mass setting, with over 140,000 DY events. Analysis of the full data 
sets has been completed and the results [71] are in good qualitative agreement with the 
high-mass data. The DY cross section ratio per nucleon for p + d to that for p + p is 
shown in Figure |^ as a function of x 2 . The acceptance of the spectrometer was largest for 
xp = x\ — x 2 > 0. In this kinematic regime the DY cross section is dominated by the 
annihilation of a beam quark with a target antiquark. To a very good approximation the 
DY cross section ratio at positive xp is given as 

vdy{p + d)/2a DY (p + p) ~ (1 + d(x 2 )/u(x 2 ))/2. (12) 

In the case that d = u, the ratio is 1. Figure |5| shows that the DY cross section per nucleon 
for p + d clearly exceeds p + p, and it indicates an excess of d with respect to u over an 
appreciable range in x 2 . 

Figure also compares the data with next-to-leading order (NLO) calculations of the 
cross section ratio using the CTEQ4M [[F|, MRS(R2) [7|, and a modified CTEQ4M parton 
distributions. The modified CTEQ4M parton distributions, in which the d + u parametriza- 
tion was maintained but d was set identical to u, were used to illustrate the cross section 
ratio expected for a symmetric d/u sea. The E866 data clearly show that d ^ u. The data 
are in reasonable agreement with the unmodified CTEQ4M and the MRS(R2) predictions 
for x 2 < 0.15. Above x 2 = 0.15 the data lie well below the CTEQ4M and the MRS(R2) 
values. 



Using an iterative procedure described in JFO , 74] , values for d/u were extracted by the 
E866 collaboration at Q 2 = 54 GeV 2 and shown in Figure [| At x < 0.15, d/u increases 
linearly with x and is in good agreement with the CTEQ4M and MRS(R2) parametrization. 
However, a distinct feature of the data, not seen in either parametrization of the parton 
distributions, is the rapid decrease towards unity of the d/u ratio beyond x 2 = 0.2. The 
result from NA51 is also shown in Figure ^] for comparison. 

The d/u ratios measured in E866, together with the CTEQ4M values for d + u, were 
used to obtain d — u over the region 0.02 < x < 0.345 (Figure [7]). As a flavor non-singlet 
quantity, d(x) — u{x) is decoupled from the effects of the gluons and is a direct measure 
of the contribution from non-perturbative processes, since perturbative processes cannot 
cause a significant d, u difference. From the results shown in Figure [!], one can obtain an 
independent determination [74] of the integral of d(x) — u(x) and compare it with the NMC 
result (Eq. [5]). E866 obtains a value Jq 1 d(x) — u(x) 
2/3 the value deduced by NMC 



— u{x) 

dx = 0.100 ± 0.007 ± 0.017, which is 



The E866 data also allow the first determination of the momentum fraction carried by 



the difference of d and u. One obtains Jqq 2 5 x d(x) — u(x) dx = 0.0065 ± 0.0010 at Q 



M is used to estimate the contributions from the unmeasured x regions, 
d(x) 



54 GeV 2 . If CTEQ< 

one finds that Jq x d(x) — u(x) dx = 0.0075 ± 0.0011. Unlike the integral of d(x) — u{x) 
the momentum integral is Q 2 dependent and decreases as Q 2 increases. The Q 2 dependence 
of the momentum fraction carried by various partons are shown in Figure The calculation 
uses both the MRS(R2) and the new MRST |67| parton distributions for comparison. 
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3.5 HERMES Semi-Inclusive Experiment 

It has been recognized for some time that semi-inclusive DIS could be used to extract 
the flavor dependence of the valence quark distributions IffSj . From quark-parton model, the 



semi-inclusive cross section, <r^, for producing a hadron on a nucleon is given by 

1 da h N (x,z) = E i eU l (x)Dl(z) 
a N (x) dz J2i^ifi( x ) 



(13) 



where D^{z) is the fragmentation function signifying the probability for a quark of flavor i 
fragmenting into a hadron h carrying a fraction z of the initial quark energy, and /$ are 
the charge and the distribution function of quark i, and <Jn{x) is the inclusive DIS cross 
section. Assuming charge symmetry for the fragmentation functions and parton distribution 
functions, one can derive the relationship 



d v {x) AR n (x) + 1 



(14) 



u v (x) A + R n (x) ' 
where 

R*(x) = (da n n + (x,z)/dz-da w n ~( y x,z)/dz)/( y daf( y x,z)/dz-daf( y x,z)/dz). (15) 
Based on a large number of semi-inclusive charged-hadron events in muon DIS from hydro- 



gen and deuterium targets, EMC extracted JFfJ the values of d v (x)/u v (x) over the range 



0.028 < x < 0.66. The EMC result agrees with neutrino measurements |77|, |78[ using a 
different method, and it demonstrates the usefulness of such semi-inclusive measurements 
for extracting information on valence quark distributions. 

Soon after the report of GSR violation by the NMC, Levelt, Mulders and Schreiber [7S| 



pointed out that semi-inclusive DIS could also be used to study the flavor dependence of sea 
quarks. In particular, 



d(x) — u(x) J(z)[l — r(x, z)] — [1 + r(x, z)} 
u(x) — d(x) J(z)[l — r(x, z)\ + [1 + r(x, z)] ' 



(16) 



where 



daf(x,z)/dz-daZ(x,z)/dz 3 1 + DZ(z)/Df(z) 

T[X ' Z) daf(x,z)/dz-daf(x,z)/dz' [Z) 5 1 - {z) / {z) ' 1 ' 

Unlike the situation for d v (x)/u v (x) which is completely independent of the fragmentation 
functions D^, Eq. |l^ shows that fragmentation functions are needed to extract the values of 
d(x) — u(x). 

Very recently, the HERMES collaboration |80| at DESY reported their measurements 



of charged hadrons produced in the scattering of a 27.5 GeV positron beam on internal 
hydrogen, deuterium, and 3 He target. The fragmentation functions DJ (z) were extracted 
from the 3 He data, while the hydrogen and deuterium data allowed a determination of 
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r(x,z). The values of (d — u)/(u — d) show no z dependence and are positive over the 
region 0.02 < x < 0.3, showing clearly an excess of d over u. Using the GRV94 LO |31 



parametrization of u(x) — d(x), the HERMES collaboration obtained d(x) — u(x) as shown 



r 0.3 
1.02 



d(x) — u(x) 



dx 



in Figure [7[ The integral of d — u over the measured x region gives J j 
0. 107 ±0. 021(stat) ±0. 017 (syst). The total integral over all x is extrapolated to be 0.16±0.03, 
consistent with the result from NMC. It is gratifying that the results from E866 and HERMES 
are in rather good agreement, even though these two experiments use very different methods 
and cover very different Q 2 values (Q 2 = 54 GeV 2 in E866 and Q 2 = 2.3 GeV 2 in HERMES). 



to situations 



It should be mentioned that semi-inclusive DIS could be extended 
involving polarized lepton beam and polarized targets in order to study the flavor dependence 
of the spin-dependent structure functions. Both the Spin Muon Collaboration (SMC) |84| 
and the HERMES Collaboration ]85| have reported the polarized valence quark distribu- 



tions, Au v (x) and Ad v (x), and the non-strange sea-quark polarization, Aq(x). 



3.6 Impact on the Parton Distribution Functions 

Since evidences for a flavor asymmetric sea were reported by the NMC and NA51, several 
groups [173, [73L performing global analysis of parton distribution function (PDF) have all 
required d to be different from u. The NMC result constrained the integral of the d — u to 
be 0.149 ± 0.039, while the NA51 result requires d/u to be 1.96 ± 0.13 at x = 0.18. Clearly, 
the x dependences of d — u and d/u were undetermined. Figures |9] and [1C] compare the 
E866 measurements of x(d — u) and d/u with the parametrizations of the MRS(R2) ||73|| , 
CTEQ4M Q7J, and GRV94 gflj, three of the most frequently used PDF's prior to E866's 



measurement. 

Recently, several PDF groups published new parametrizations taking into account of new 
experimental information including the E866 data. The parametrization of the x dependence 
of d—u is now strongly constrained by the E866 and HERMES data. In particular, d(x)—u(x) 
or d(x)/u(x) are parametrized as follows; 

MRST 



d{x) - u{x) = 1.29x ai83 (l - x) 9 ' 808 (l + 9.987x - 33.34x 2 ), Q 2 = 1 GeV 2 , 
GRV98 @: 

d(x) - u(x) = 0.20a;- ' 57 (l - a;) 12 - 4 (l - 13.3x a5 + QO.Ox), Q 2 = 0.4 GeV 2 , 
CTEQ5M [|7|, ||: 

d(x)/u(x) = 1 - 1.095x + 3.159x a5 /[l + {(x - 0.188)/0.116) 2 ], Q 2 = 1 GeV 2 . (18) 

As shown in Figure [I~0|, these new parametrizations give significantly different shape for d/u 
at x > 0.15 compared to previous parametrizations. Table [j] also lists the values of d — u 
integral from various recent PDF's. 

It is interesting to note that the E866 data also affect the parametrization of the valence- 
quark distributions. Figure |TT] shows the NMC data for F% — F£ at Q 2 = 4 GeV 2 , together 
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with the fits of MRS(R2) and MRST. It is instructive to decompose F$ (x) — F£{x) into 
contributions from valence and sea quarks: 

1 2 r_ - i 

F$ (x) — F%(x) = -x [u v (x) — d v (x)] + -x u{x) — d(x) . (19) 



As shown in Figure |TTJ, the E866 data provide a direct determination of the sea-quark 
contribution to F^ — F^- In order to preserve the fit to F^ — F^, the MRST's parametrization 
for the valence-quark distributions, u v — d v , is significantly lowered in the region x > 0.01. 
Indeed, one of the major new features of MRST is that d v is now significantly higher than 
before at x > 0.01. Although the authors of MRST attribute this to the new W^-asymmetry 
data from CDF |8{| and the new NMC results on F^/F^ @, it appears that the new 
information on d(x) — u(x) has a direct impact on the valence-quark distributions too. 

Another implication of the E866 data is on the behavior of Ff — F 2 n at small x. In 
order to satisfy the constraint J [u v (x) — d v (x)]dx = 1, the MRST values of u v (x) — d v (x) 
at x < 0.01 are now much larger than in MRS(R2), since u v (x) — d v (x) at x > 0.01 are 
smaller than before. As a consequence, Ff — Fg is increased at small x and MRST predicts 
a large contribution to the Gottfried integral from the small-x {x < 0.004) region, as shown 
in Figure [T2. If the MRST parametrization for Ff — F£ at x < 0.004 were used, NMC would 
have deduced a value of 0.252 for the Gottfried integral, which would imply a value of 0.122 
for the d — u integral. This would bring better agreement between the E866 and the NMC 
results on the d — u integral. 



3.7 Comparison Between Various Measurements 

Are the various measurements of sea quark flavor asymmetry consistent among them? 
In particular, is the E866 result consistent with the earlier E772 and the NA51 DY experi- 
ments, and with the HERMES and NMC DIS measurements? To address this question, we 
first compare E866 with E772, both of which are DY experiments using 800 GeV proton 
beam with essentially the same spectrometer. Figures || and [| show that the MRST parton 
distributions, which determined d — u based on the E866 data, can also describe the E772 
data very well, and we conclude that the E866 and E772 results are consistent. 

Although both NA51 and E866 measured a(p + d) /2a(p + p) to extract the values of d/u, 
some notable differences exist. As mentioned earlier, NA51 measured the ratio at a single 
value of X2 (22 = 0.18) near xp ~ using 450 GeV proton beam, while E866 used 800 GeV 
proton beam to cover a broader range of X2 at xp > 0. It is instructive to compare the NA51 
result at X2 = 0.18 with the E866 data at X2 = 0.182. Table [| lists the kinematic variables 
and physics quantities derived from these two data points. It is interesting to note that the 
values of a(p + d)/2a(p + p) at X2 = 0.18 are actually very similar for NA51 and E866, even 
though the derived values for d/u differ significantly. This reflects the difference in for 
both experiments, making the values of d/u extracted from a(p + d) /2a(p + p) different. The 
other difference is Q 2 , being a factor of 3.6 higher for E866. Using MRST and CTEQ5M to 
estimate the Q 2 dependence of d/u, we find that the NA51 value of d/u is reduced by ~ 3 
% going from Q 2 = 27.2 GeV 2 to Q 2 = 98.0 GeV 2 . This brings slightly better agreement 
bewteen NA51 and E866. 
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The methods used by HERMES and E866 to determine d — u are totally different, and 
it is reassuring that the results came out alike, as shown in Figure |7|. The d — u values from 
HERMES are in general somewhat larger than those of E866. At a relatively low mean Q 2 
of 2.3 GeV 2 , the HERMES experiment could be subject to high-twist effect |90fl . Additional 
data from HERMES are expected to improve the statistical accuracy. 

The comparison between E866 and NMC in terms of the integral of d — u has been dis- 
cussed earlier. A possible origin for the apparent differences of the integral was also discussed 
in Section 3.6. 

4. ORIGINS OF THE d/u ASYMMETRY 
4.1 Pauli Blocking 

The earliest experiments indicated that the value of the Gottfried integral might be less 
than 1/3, leading to speculation regarding the origin of this reduction. Field and Feynman 
suggested |24| that it could be due to Pauli blocking in so far as uu pairs would be suppressed 
relative to dd pairs because of the presence of two u-quarks in proton as compared to a single 
(i-quark. Ross and Sachrajda [|5f| questioned that this effect would be appreciable because of 
the large phase-space available to the created qq pairs. They also showed that perturbative 
QCD would not produce a d, u asymmetry. Steffens and Thomas |91j recently looked into 



this issue, explicitly examining the consequences of Pauli blocking. They similarly concluded 
that the blocking effects were small, particularly when the antiquark is in a virtual meson. 

The small d, u mass difference (actually, > m u ) of 2 to 4 MeV compared to the 
nucleon confinement scale of 200 MeV/c does not permit any appreciable difference in their 
relative production by gluons. At any rate, one observes a surplus of d which is the heavier 
of the two species. As pointed out above, blocking effects arising from the Pauli exclusion 
principle should also have little effect. Thus another, presumably non-perturbative, mecha- 
nism must be found to account for the large measured d, u asymmetry. 

4.2 Meson-Cloud Models 

A natural origin for this flavor asymmetry is the virtual states of the proton containing 
isovector mesons. This point appears to have first been made by Thomas in a publication 



92] treating SU(3) symmetry breaking in the nucleon sea. Sullivan JS^j had earlier shown 
that virtual meson-baryon states directly contribute to the nucleon's structure function. A 
large number of authors have contributed to calculating the asymmetry from this perspective, 



so recent reviews |)4|, ^] should be consulted for a complete list of contributions. 

Conservation of electric charge and isospin naturally creates a flavor asymmetry when 
isovector mesons are involved. For example the Fock decomposition of the proton into its 
ttN components yields 

\P)^\[1 bo7r°) + y||n 7r + }, (20) 
where po and are regarded as proton and neutron states with symmetric seas. The d/u 



16 



ratio for the configurations shown in Eq. is 5! Thus the presence of such virtual states can 
readily generate a large flavor asymmetry. Allowing for the possibility of it A configurations 
it is easy to show that the integrated d, u difference in the proton is given by 

f \d(x) -u(x)]dx = -(2a -b), (21) 
Jo 3 

where a(b) is the probability of 7riV(7rA) component in the proton. The value extracted from 
E866 for the above integral is 0.100 ±0.018, leading to a = 0.2 ±0.036 if one accepts b ~ a/2 
as found in many calculations. 



Many attempts |JJ| [J7|, [51], |10U| , |101| , |102j | have been made to calculate the flavor 



asymmetry due to isovector mesons. Most start with the following convolution expressions: 

\p) = a \ Po ) + J2a P MB\M)\B), (22) 



MB 



xq p (x, Q 2 ) = J2 a liB / dy fMB(y) - Qm(-,Q 2 ), (23) 

mb y y 



where 



r . . 9m p b [ tmin F(t,m p ,m B ) 2 2 i u ,,,,, 
JmbVU) = / dt ~u T^r F M P B\t, A )> tmin = m p y - m B - . (24) 

1D7T J —oo [Z — TYl^ J 1 — y 



In the above expressions x is the fraction of proton's momentum carried by the antiquark, and 
y is fraction carried by the meson (M). The meson-proton-baryon couplings are characterized 
by coupling constants gM P B, and form factors F MpB (t,A) where A is a cutoff parameter. 
F(t, m p , ms) is a kinematic factor depending on whether B is in the baryon octet or decuplet. 
As pions are the only mesons usually considered and the baryons are usually restricted to 
nucleons and deltas, the coupling constants are well known and the partonic structure of the 
pion, q n (x,Q 2 ), is fixed by measurement of the DY process using high energy pion beams. 
The only uncertainties are the form factors F wpN {t) and F np &(t). One attempts to determine 
these form factors by using | |103| , |104| | the measured yields from a variety of high energy 



hadronic reactions at small pt such as p + p — > n + X and p ± p — ► A ++ ± X. Even though 
there is a sizable amount of available data, employing such a procedure does not produce a 
precise result. The cutoff parameters used in the extracted form factors are of the order of 1 
GeV but the uncertainties in their values produce factors of two in the predicted antiquark 
content of the nucleon. 

Even though calculated value of the integral of d(x) — u(x) is often in agreement with 
experiment, it is more difficult to achieve a quantitative fit to the measured x dependence 
of the difference [74]. Figure |7| compares d(x) — u(x) from E866 with a pion-cloud model 



calculation, following the procedure detailed by Kumano []97[. A dipole form, with A = 1.0 



GeV for the nNN form factor and A = 0.8 GeV for the 7rA^A form factor, was used. 
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Calculations of d(x)/u(x) are even more unsuccessful, as knowledge of the x dependence of 
the symmetric sea is required in this instance JT4[ . 



It is instructive to compare the pion-model prediction with the current PDF parametriza- 
tion of x(d + u). Figure [13] shows that at small x (x < 0.1) the valence quarks in the pion 
cloud account for less than 1/3 of the d + u content in the proton. In contrast, at large x 
(x > 0.5) the pion model would attribute all of d + u to the pion cloud. 

More recent attempts 



ioa, 104 



to calculate the asymmetry due to isovector mesons 
find a smaller 7rA component in the nucleon than is presented following Eq. |2T| The ttN 
component is still about 20% while the irA piece is down around 6%. As a result the inte- 
grated asymmetry is too large, typically 0.16. However considerable progress has been made 
by Reggeizing the virtual mesons [|104j| . This procedure has two interesting consequences, 
first it shows that the principal contributions to the asymmetry come from virtual pions 
and secondly it reduces the contributions from these pions at large x, allowing the ratio 
d(x)/u(x) to approach 1 for x > 0.3. Another recent work |105| describes meson-baryon 
fluctuation through gluon splitting and a phenomenological recombination mechanism, and 
it also reproduces the x dependence of d/u well. 

Recently, the flavor asymmetry of the nucelon sea was computed in the large- N c limit, 
where the nucleon is described as a soliton of an effective chiral theory ||106| , |107|| . In this 
chiral quark-soliton model, the flavor non-singlet distribution, u(x) — d(x), appears in the 
next-to-leading order of the 1/N C expansion ||108| , |109|| . The E866 d(x) — u(x) data were 



shown to be well described by this model 107 



4.3 Chiral Models 



An alternative approach [^B, |110| , |111| , |112| , |113| , |114| , |115|| also employing virtual pions 



to produce the d, u asymmetry uses constituent quarks and pions as the relevant degrees of 
freedom. Such models are usually referred to as chiral models [ITB|. In this model, a portion 



of the sea comes from the couplings of Goldstone bosons to the constituent quarks, such as 
u — > dn + and d — > utt~. The excess of d over u is then simply due to the additional up 
valence quark in the proton. 

The chiral models have less dynamical freedom, and the meson-baryon expansion can 
be shown to predict d(x)/u(x) < 11/7 independent of x or Q 2 . This limit is exceeded for 
0.12 < x < 0.2 as shown in Figure |[ The limit on the maximum value of the ratio can 
be traced to the fact that there is no mechanism to suppress virtual 7rA-like configurations 
which are seen in Eq. [H] to reduce the asymmetry. Some extensions of the simple chiral 
quark model can obtain agreement with certain features of the observed asymmetry but 



only with the introduction of additional parameters. A calculation |74 , |113|| of d(x) — u(x) 
(Figure |7|) yielded too soft a distribution relative to experiment, indicating |73| that more 
dynamics must be included in the chiral quark model if it is to produce more than rough 
qualitative agreement with the observed asymmetry in the up, down sea of the nucleon. 



4.4 Instanton Models 
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Instantons have been known as theoretical constructs since the seventies [ IT7 , 118 , 119 |. 



They represent non-perturbative fluctuations of the gauge fields that induce transitions be- 
tween degenerate ground states of different topology. In the case of QCD, the collision 
between a quark and an instanton flips the helicity of the quark while creating a qq pair 
of different flavor. Thus, interaction between a u quark and an instanton results in a u 
quark of opposite helicity and either a dd or ss pair. Such a model has the possibility of 
accounting for both the flavor asymmetry and the "spin crisis" ||120| , |121|| . However, the 
model has proven difficult to exploit for this purpose. There is only one case [ 122 | of its 
being employed to explain these anomalous effects. In the case of the d, u flavor asymmetry 
the authors of ref. [ |122j | fit the instanton parameters to reproduce the violation of the GSR 
observed by NMC. The prediction [122] at large x, d(x)/u(x) — > 4, is grossly violated by 
experiment (see Figure |^). Thus, it appears that while instantons have the possibility for 
accounting for flavor and spin anomalies, the approach is not yet sufficiently developed for a 
direct comparison. The final state created via an instanton collision is quite similar to that 
created via the emission of meson in the chiral model. 



4.5 Charge Symmetry Violation 

Charge symmetry is believed to be well respected in strong interaction. Extensive exper- 
imental searches for charge symmetry violation (CSV) effects in various nuclear processes 
reveal an amount on the order of 0.3% ||123|| . This is consistent with the expectation that CSV 
effects are caused by electromagnetic interaction and by the small mass difference between 
the u and d quarks 



It has been generally assumed that the parton distributions in hadrons obey charge 
symmetry. This assumption enables one to relate the parton distributions in the proton and 
neutron in a simple fashion, u p (x) = d n (x),d p (x) = u n (x), etc. Indeed, charge symmetry is 
usually assumed in the analysis of DIS and DY experiments, which often use nuclear targets 
containing both protons and neutrons. Charge symmetry is also implicit in the derivation of 
many QCD sum rules, including the Gottfried sum rule, the Adler sum rule, and the Bjorken 
sum rule. 

The possibility that charge symmetry could be significantly violated at the parton level 
has been discussed recently by several authors JT25], [126], [12?|, [128], [129], [130], [13T], [132|. Ma 
and collaborators [ |125| , |126|] pointed out that the violation of the GSR can be caused by CSV 



as well as by flavor asymmetry of the nucleon sea. They also showed that DY experiments, 
such as NA51 and E866, are subject to both flavor asymmetry and CSV effects. Forte 
estimated | |128| | the relative size of the CSV and flavor asymmetry via a combined analysis 
of the GSR violation and the nucleon a term. Using the bag model, Rodionov et al [|130| 
showed that a significant CSV effect of ~ 5% could exist for the "minority valence quarks" 
[i.e. d p (x) and u n (x)} at large x (x > 0.4). A model study [|132|| of CSV for sea quarks shows 
that the effect is very small, roughly a factor of 10 less than for valence quarks. Londergan 
& Thomas have reviewed the role of CSV for parton distributions ||133|| . 

Recently Boros, Londergan and Thomas (BLT) ||134j| compared the parton distributions, 
x,Q 2 ) and F^^x^Q 2 ). The F^^x^Q 2 ) structure function extracted from DIS muon 



F uN 
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scattering is defined as 



Ff(XtQ ,^mam±m^ (25) 



5 2 8 

—x[u + u + d + d+-(s + s) + -(c + c)]. (26) 



In going from Eq. to EB, the (x, Q ) dependence has been suppressed and charge symmetry 



(w p = d n = u, etc.) has been employed. For neutrino DIS one has 

F% N (x, Q 2 ) =x[u + u + d + d + 2s + 2c], (27) 

again using charge symmetry and suppressing the (x, Q 2 ) dependence in the parton distri- 
butions. Neglecting the contribution of charmed quarks in the nucleon and correcting for 
the small difference due to strange quark contributions one expects 

^F£ N (x,Q 2 ) , 

^FW K L (28) 

After making the necessary small corrections, BLT found that the ratio in Eq. ^ at 
common values of (x, Q 2 ) is 1 for x > 0.1 but for x < 0.1 the ratio progressively decreases 
below 1 as x decreases. They suggested that this effect might be due to a violation of charge 
symmetry in the sea of the nucleon. To achieve agreement with experiment, BLT found it 
necessary to set d n (x) ~ 1.25u p (x) and u n (x) ~ 0.75d p (x) at small x. As the total number 
of sea quarks is kept equal in the neutron and proton (u p (x) + d p (x) ~ u n (x) + d n (x)), one 
deduces d p (x) ~ u p (x) and d n (x) ~ (5/3)u n (x), a most striking asymmetry between the 
proton and the neutron sea. 

How would this claim of large CSV affect the E866 analysis of the flavor asymmetry? 
First, CSV alone could not account for the E866 data. In fact, an even larger amount of 



flavor asymmetry is required to compensate for the possible CSV effect ||134|| . Second, there 
has been no indication of CSV for x > 0.1. Thus, the large d/u asymmetry from E866 for 
x > 0.1 is not affected. Fortunately, this radical alteration of the conventional view does 
not appear to be the CclSC, clS shown recently by Bodek et al. [ |135|| , who showed that the 
measured asymmetry of iy + (iy~) production at CDF is consistent with charge symmetry 
and in strong disagreement with the suggestion of BLT. The reasons for the discrepancy of 



the DIS data with Eq. 28 must be found elsewhere 



5. FURTHER IMPLICATIONS OF THE MESON-CLOUD MODELS 

5.1 Strange Sea of the Nucleon 

Models in which virtual mesons are admitted as degrees of freedom have implications that 
extend beyond the d, u flavor asymmetry addressed above. They create hidden strangeness 
in the nucleon via such virtual processes as p — > A + K + , E + K, etc. Such processes are of 
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considerable interest as they imply different s and s parton distributions in the nucleon, a 
feature not found in gluonic production of ss pairs. This subject has received a fair amount 
of attention in the literature ||103| , |136| , p. 37] , |138| , |139|| but experiments have yet to clearly 
identify such a difference. Thus in contrast to the d, u flavor asymmetry, to date there is no 
positive experimental evidence for ss contributions to the nucleon from virtual meson-baryon 
states PLIHOU- 

A difference between the s and s distribution can be made manifest by direct measurement 
of the s and s parton distribution functions in DIS neutrino scattering, or in the measurement 
of the q 2 dependence of the strange quark contribution (-Ff s (g 2 )) to the proton charge form 
factor. This latter case is not well known and follows from a suggestion of Kaplan and 
Manohar [|141|| regarding the new information contained in the weak neutral current form 



factors of the nucleon. Measurement of these form factors allows extraction of the strangeness 
contribution to the nucleon's charge and magnetic moment and axial form factors. The 
portion of the charge form factor Fl s (q 2 ) due to strangeness clearly is zero at q 2 = 0, but if 
the s and s distributions are different the form factor becomes non-zero at finite q 2 . These 
"strange" form factors can be measured in neutrino elastic scattering ||142|| from the nucleon, 
or by selecting the parity-violating component of electron-nucleon elastic scattering, as is 
now being done at the Bates ||143|| and Jefferson Laboratories [|144j . 

It is worth pointing out that there is a relationship between the parton distributions and 
the form factors of a hadron. If the neutron's charge form factor is explained in terms of a 
particular meson-baryon expansion, then one should expect that the expansion is consistant 
with the neutron's partonic structure. Little work appears to have been done bringing these 
descriptions together. Below is an example showing the impact of the nucleons pionic con- 
tent, inferred from the flavor asymmetry, on the spin distribution in the nucleon. 



5.2 Spin Dependent Structure Functions 

As the pion is a J 71 " = 0~ hadron it must be emitted as a p-wave requiring that the 
nucleon flip its spin 1/3 of the time upon emitting a pion. Thus, it appears that such a 
process might readily account for the reduction of g A from 1.667 and the overall reduction 
in the quark contribution to nucleon spin as observed in DIS. The contribution to g A from 
virtual pions is 

5 20 32 

g A = Ait - Ad = - - -{2a + b) + -V2ab (29) 

= 1.53 for a = 0.2,6 = 0.1, (30) 

and to the total nucleon spin due to up and down quarks 

2 

Aq = Au + Ad = 1 - -(2a - 6) (31) 

3 

= 0.80 for a = 0.2, b = 0.1. (32) 

Thus, virtual pions do reduce the fraction of nucleon spin carried by up and down quarks but 
the effect is not sufficient to reduce g A from its SU(6) value of 5/3 to the experimental value of 
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1.256. Eqs. and [3D], together with the values of d—u integral determined from experiments, 
allow us to calculate g& as a function of a. The different values of the d — u integral from the 
NMC and E866 correspond to two different curves shown in Figure [T3|. Figure |T^ illustrates 



the possible range of g& achievable from virtual pions alone. Presumably relativistic effects 
are responsible for the bulk of that reduction |145| . Similarly, virtual pions reduce the total 



spin carried by u and d quarks but not enough to be a dominant effect as Aw + Ad « 0.25. 

There is an aspect of spin distributions predicted by virtual pion emission that is present 
in the data. The spin carried by antiquarks in the virtual pion model must be zero. Indeed 



the measured value f84| , Aq = 0.01 ± 0.04 ± 0.03 is consistent with that picture. It is very 
likely that meson will prove to play an important part in understanding the nucleon's spin 
structure but that role is not clear at present. 

5.3 Sea Quark Distributions in Hyperons 

Dilepton production using meson or hyperon beams offers a means of determining par- 
ton distributions of these unstable hadrons. Many important features of nucleon parton 
distributions, such as the flavor structure and the nature of the non-perturbative sea, find 
their counterparts in mesons and hyperons. Information about meson and hyperon parton 
structure could provide valuable new insight into nucleon parton distributions. Furthermore, 
certain aspects of the nucleon structure, such as the strange quark content of the nucleon, 
could be probed with kaon beams. 

No data exist for hyperon- induced dilepton production. The observation of a large d/u 
asymmetry in the proton has motivated Alberg et al. ||146| , |147|| to consider the sea-quark 



distributions in the E. The meson-cloud model implies a d/u asymmetry in the S + even 
larger than that of the proton. However, the opposite effect is expected from SU(3) sym- 
metry. Although relatively intense S + beams have been produced for recent experiments at 
Fermilab, this experiment appears to be very challenging because of large pion, kaon, and 
proton contaminations in the beam. 

6. FUTURE PROSPECTS 

6.1 d/u at Large and Small x 

The interplay between the perturbative and non-perturbative components of the nucleon 
sea remains to be better determined. Since the perturbative process gives a symmetric d/u 
while a non-perturbative process is needed to generate an asymmetric d/u sea, the relative 
importance of these two components is directly reflected in the d/u ratios. Thus, it would 
be very important to extend the DY measurements to kinematic regimes beyond the current 
limits. 

The new 120 GeV Fermilab Main Injector (FMI) and the proposed 50 GeV Japanese 
Hadron Facility ||148|| (JHF) present opportunities for extending the d/u measurement to 
larger x (x > 0.25). For given values of x\ and X2 the DY cross section is proportional to 
1/s, hence the DY cross section at 50 GeV is roughly 16 times greater than that at 800 



GeV! Figure [15] shows the expected statistical accuracy for a(p + d)/2a(p + p) at JHF ||149 
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compared with the data from E866 and a proposed measurement [ 15U|| using the 120 GeV 



proton beam at the FMI. A definitive measurement of the d/u over the region 0.25 < x < 0.7 
could indeed be obtained at FMI and JHF. 

At the other end of the energy scale, RHIC will operate soon in the range 50 < y/s < 500 
GeV/nucleon. The capability of accelerating and colliding a variety of beams from p + p, 
p + A, to A + A at RHIC offers a unique opportunity to extend the DY d/u measurement 
to very small x. 

6.2 W Production 

To disentangle the d/u asymmetry from the possible CSV effect, one could consider W 
boson production, a generalized DY process, in p + p collision at RHIC. An interesting 
quantity to be measured is the ratio of the p + p — > W + + X and p + p — > W~ + X cross 



sections |[L51|| . It can be shown that this ratio is very sensitive to d/u. An important fea- 
ture of the W production asymmetry in p + p collision is that it is completely free from 
the assumption of charge symmetry. Figure [16| shows the predictions for p + p collision at 



y/s = 500 GeV. The dashed curve corresponds to the d/u symmetric MRS SO' ||152f| structure 
functions, while the solid and dotted curves are for the d/u asymmetric structure function 
MRST and MRS(R2), respectively. Figure [16] clearly shows that W asymmetry measure- 
ments at RHIC could provide an independent determination of d/u. 

6.3 Strange Sea in the Nucleon 

As discussed earlier, an interesting consequence of the meson-cloud model is that the s and 
s distributions in the proton could have very different shapes, even though the net amount 
of strangeness in the proton vanishes. By comparing the v and v induced charm production, 
the CCFR collaboration found no difference between the s and s distributions |140fl . More 



precise future measurements would be very helpful. Dimuon production experiments using 
K ± beams might provide an independent determination of the s/s ratio of the proton, 
provided that our current knowledge on valence-quark distributions in kaons is improved. As 
discussed in Section 5.1, ongoing measurements of Ff s via parity- violating electron-nucleon 
scattering should shed much light on the possible difference between s and s distributions. 

6.4 Sea Quark Polarization 

Polarized DY and W ± production in polarized p + p collision are planned at RHIC ||153| 



and they have great potential for providing qualitatively new information about antiquark 
polarization. At large xp region [xp > 0.2), the longitudinal spin asymmetry All i n the 
p + p DY process is given by | ]154| , |155| | 

Aw 

A%(x u x 2 ) w g x {x 1 )/F 1 {x 1 ) x — (x 2 ), (33) 

u 

where gi(x) is the proton polarized structure function measured in DIS, and Au(x) is the 
polarized u distribution function. 

Eq. shows that u polarization can be determined using polarized DY at RHIC. Addi- 



23 



tional information on the sea-quark polarization can be obtained via W production ||156 



The parity-violating nature of W production implies that only one of the two beams need 



to be polarized. At positive xf (along the direction of the polarized beam), one finds[ |156|j , 



AY+m—M, and 4T-^W (34) 

where is the single-spin asymmetry for W production. Eq. |34| shows that the flavor 
dependence of the sea-quark polarization can be revealed via W ± production at RHIC. A 
remarkable prediction of the chiral quark-soliton model is that the flavor asymmetry of po- 
larized sea-quark, Au(x) — Ad(x), is very large ||157|| . This is in striking contrast to the 
meson cloud model which predicts very small values for Au(x) — Ad(x) [ 158 , T5§ \. Future 



DY and W ± production experiments at RHIC could clearly test these models ||160 



7. CONCLUSION 

The flavor asymmetry of the nucleon sea has been clearly established by recent DIS and 
DY experiments. The surprisingly large asymmetry between u and d is unexplained by per- 
turbative QCD. Thus far, three distinct explanations have been offered in the literature for 
the origin of the observed d, u asymmetry. The first is Pauli blocking which, while appealing, 
is difficult to calculate and appears to produce too small an effect to be the sole origin of the 
large observed asymmetry. The second involves virtual isovector mesons, mostly pions, in the 
nucleon. Such descriptions necessarily require non-perturbative QCD and, apart from lat- 
tice gauge calculations, demand additional parameters and possess systematic uncertainties. 
However, as these virtual mesons readily generate a large d, u asymmetry, many authors, 
using a variety of techniques to invoke and justify their approachs, have obtained qualitative 
agreement with the experimental measurements. The third explanation involves instantons, 
but the theory is not sufficiently developed to allow quantitative comparison to the asymme- 
try data. Future experiments will test and refine these models. They will further illuminate 
the interplay between the perturbative and non-perturbative nature of the nucleon sea. 
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Table 1: Values of the integral ^{d{x) — u{x))dx for various parton distribution functions. 



PDF 


integral 


CTEQ4M 


0.108 


MRS(R2) 


0.162 


GRV94 


0.163 


CTEQ5M 


0.124 


MRST 


0.102 


GRV98 


0.126 
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Table 2: Comparison between the 450 GeV/c NA51 result and the 800 GeV/c E866 data 
point near x 2 = 0.18. 





(*2> 


(x F ) 


(M^XGeV) 




d/u 


NA51 


0.18 


0.0 


5.2 


1.099 ± 0.039 


1.96 ± 0.246 


E866 


0.182 


0.192 


9.9 


1.091 ± 0.044 


1.41 ± 0.146 



0.12 



0. 



0.08 



0.06 



0.04 - 



0.02 




Figure 1: EMC [f|| and SLAC measurements of F 2 P 
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Figure 2: NMC measurements of F% — F£ (solid data points) and the Gottfried integral 



(open data points). The triangular data points correspond to results published in 1991 ||6 



while the circular data points represent a more recent analysis in 1994 |5{J. The extrapolated 
value of Gottfried integral (Sq) and the expected GSR value are also indicated. 
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Figure 3: The ratios of W over isoscalar targets DY cross sections from E772 [64|. The 



dashed curve is a calculation using the d/u asymmetric parton distributions suggested in 



Ref. [65). The solid curve corresponds to a calculation using the recent MRST distribution 
functions. 
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Figure 4: The p+ 2 H DY differential cross sections from E772 ||64|| . The dotted (dashed) curve 
is a calculation using the parton distributions from Ref. |)S with (without) d/u asymmetry. 
The solid curve uses the MRST distribution functions. These are leading-order calculations 
normalized to the data. 
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Figure 5: The ratio a pd /2a pp of Drell-Yan cross sections vs. x 2 . The curves are next-to- 
leading order calculations, weighted by acceptance, of the Drell-Yan cross section ratio using 
the CTEQ4M and MRS(R2) parton distributions. Also shown is a leading-order calculation 
using CTEQ4M (dotted). In the lower CTEQ4M curve d — u has been arbitrarily set to as 
described in the text. The errors are statistical only. There is an additional 1% systematic 
uncertainty common to all points. 
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Figure 6: The ratio of d/u in the proton as a function of £2 extracted from the Fermilab 
E866 cross section ratio. The curves are from various parton distributions. The error bars 
indicate statistical errors only. An additional systematic uncertainty of ±0.032 is not shown. 
Also shown is the result from NA51, plotted as an open box. 
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Figure 7: Comparison of the E866 [70| d — u results at Q 2 = 54 GeV 2 /c 2 with the predictions 
of pion-cloud and chiral models as described in the text. The data from HERMES |8(J are 
also shown. 




Figure 8: Q-dependence of the proton's momentum carried by various partons calculated 
using the MRS(R2) and MRST parton distribution functions. 




Figure 9: Comparison between the x(d — u) results from E866 with the parametrizations of 
various parton distribution functions. The dashed (solid) curves correspond to PDFs before 
(after) the E866 results were obtained. 
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Figure 10: Comparison between the d/u results from E866 with the parametrizations of 
various parton distribution functions. The dashed (solid) curves correspond to PDFs before 
(after) the E866 results were obtained. 
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Figure 11: F$ — F£ as measured by NMC at Q 2 = 4 GeV 2 compared with predictions based 
on the MRS(R2) and MRST parametrizations. Also shown are the E866 results, evolved 
to Q 2 = 4 GeV 2 , for the sea-quark contribution to F$ — F£. For each prediction, the top 
(bottom) curve is the valence (sea) contribution and the middle curve is the sum of the two. 
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Figure 12: F$ — F£ as measured by NMC at Q 2 = 4 GeV 2 compared with parametrization of 
MRST. The dashed curve corresponds to 0.21x 62 , a parametrization assumed by the NMC 
for the unmeasured small-x region when the Gottfried integral was evaluated. 
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Figure 13: Comaprison of x(d + u) obtained from the valence quarks in the pion cloud with 
the parametrization of the MRST parton distribution functions. 
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Figure 14: as a function of a determined from Eqs. [21] and [30|. The upper and lower solid 
curve is obtained by using the values d — u determined by E866 and NMC, respectively. The 
values of gA for the SU(6) limit and from the experiments are also indicated by the dashed 
curves. 



47 
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Figure 15: Projected statistical accuracy for a^j2a v in a 100-day run at JHF. The E866 data 
and the projected sensitivity for a proposed measurement ||150|| at the 120 GeV Fermilab 
Main-Injector are also shown. 
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Figure 16: Predictions of a(p + p — > W + X)/a(p + p — > W~X) as a function of x F at = 
500 GeV. The dashed curve corresponds to the d/u symmetric MRS SO' structure functions, 
while the solid and dotted curves are for the d/u asymmetric structure function MRST and 
MRS(R2), respectively 



